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ABSTRACT 

The recent spectacular progress in the experimental and theoretical understanding 
of graphene, the basic constituent of graphite, is applied here to compute, from first 
principles, the UV extinction of nano-particles made of stacks of graphene layers. The 
theory also covers cases where graphene is affected by structural, chemical or orienta- 
tion disorder, each disorder type being quantitatively defined by a single parameter. 
The extinction bumps carried by such model materials are found to have positions 
and widths falling in the same range as the known astronomical 2175 A features: as 
the disorder parameter increases, the bump width increases from 0.85 to 2.5 /im~^, 
while its peak position shifts from 4.65 to 4.75 /im~^. Moderate degrees of disorder 
are enough to cover the range of widths of the vast majority of observed bumps (0.75 
to 1.3 /im~^). Higher degrees account for outliers, also observed in the sky. 

The introduction of structural or chemical disorder amounts to changing the initial 
sp'^ bondings into sp^ or sp^ , so the optical properties of the model material become 
similar to those of the more or less amorphous carbon-rich materials studied in the 
laboratory: a-C, a-C:H, HAC, ACH, coals etc. The present treatment thus bridges 
gaps between physically different model materials. 

Key words: astrochemistry — ISM: lines and bands — dust 



1 INTRODUCTION 

The ISEC (Interstellar Extinction Curve) is dominated, in 
the vis/UV, by a strong and ubiquitous feature, peaking at 
5.7 eV (or 4.6 /xm~^). The discoverers of this feature sug- 
gested that it might be carried by small crystalline graphite 
particles (Stecher and Donn 1965). Pure bulk crystalline 
graphite in the laboratory displays an extinction feature 
near 4.25 eV, but, in fine powder form, its extinction peaks 
at 5.7 eV and was assigned to a "surface plasmon mode", 
or Frohlich resonance, due to the tv electron of graphite 
(Bohren and Huffman, 1983); it is also called "bump" in 
astronomer parlance. Since then, however, together with de- 
velopments in this direction, several other candidate carriers 
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have been predicated, ranging from molecular to solid, more 
or less hydrogenated and more or less graphitic. A necessar- 
ily very short and partly subjective sample of the latest of 
these models includes mixtures of PAHs (Cecchi-Pestellini et 
al. 2008); Bucky Onions (fullerenes) ; fullerenes/buckyonions 
models (e.g. De Heer and Ugarte 1993, Wada and Toku- 
naga 2006, Chhowalla et al. 2003, Ruiz et al. 2005, Li et al. 
2008); mixtures of stacks of PAH molecules and HAC/a-C:H 
(hydrogenated amorphous carbon (Duley and Seahra 1999); 
mixtures of ellipsoidal amorphous carbon grains (Mennella 
et al. 1998); agglomerated HAC grains (Schnaiter et al. 
1998 ), and PCG (polycrystalline graphite, Papoular and 
Papoular 2009). The common purpose is to satisfy both 
tight constraints imposed by astronomical observations: rel- 
atively large variations of width (roughly 0.7 to 1.3 /im~^) 
while the peak position remains fixed at 4.6 /xm"^ to ^2 
% (see Fitzpatrick and Massa 1986, 2007). Figure 16 of 
Fitzpatrick and Massa p007p plots all observed features in 
a graph of the peak wavenumber, z/q, against the FWHM 
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(Full Width Half- Maximum) width, 7. It shows that the 
vast majority of representative points lie within a clus- 
ter centered at {^-^ 0.9 /xm""^, 4.59 /xm""^), with dimensions 
(^ 0.3 /im~^, 0.1 /im~^). However, many outliers are also 
observed, which seem to extend this cloud towards larger 
widths, with rapidly decreasing density, as confirmed by the 
histograms of fig. 17 of the same reference (see also, for in- 
stance, Cardelli and Savage 1988). 

A large number of carbonaceous materials are found to 
exhibit the tt and a UV resonances of graphite with various 
relative intensities, even if highly loaded with hydrogen. In 
general, the former decrease considerably as structural dis- 
order and chemical impurities increase. This is illustrated 
by artificially produced a-C:H (see for instance Fink et al. 
1984) as well as natural coals extracted from mines of de- 
creasing depth, i.e. of decreasing degree of graphitization. 
Papoular et al. (|1993p and Papoular et al. (1995) showed 
that the tt feature is exhibited not only by graphite, but also 
by coals with decreasing degree of graphitization, until it is 
completely damped by excessive H content and sp^ bond- 
ings. The inverse evolution can be induced, starting from 
the less graphitized coal or kerogen, by thermal annealing 
and/or high energy irradiation. Mennella et al. (|1995p and 
Mennella et al. (|1998p , studying artificially produced car- 
bonaceous materials, also observed and used this general 
property of carbon based materials. 

For all types of materials, it was found that the a fea- 
ture can still be observed when the tt feature has long been 
drowned. This is to be expected since the former is associ- 
ated with ppa types of bonding (in which principal axes of 
coupled p orbit als lie along the same line; see Castro Neto et. 
al. (^2009J), which survive the transition to sp^ hybridization 
that dominates structurally disordered materials. Robertson 
(|1986p , in his study of HAC, compared extensively the op- 
tical and other properties of graphite with those of disor- 
dered carbons (HAC, glassy carbon, a-C, a-C:H, etc.). It is 
striking to see the kind of continuity linking these various 
carbon-rich materials and their transformation from one into 
another through irradiation and/or thermal annealing. Sev- 
eral attempts at modelling these disordered materials were 
made, with a view to understanding the physics underlying 
their spectral properties. For instance, Stenhouse and Grout 
(|1978p proposed the "domain model" , which pictures a-C as 
microcrystallites of graphite inter-linked by an sp^-bonded 
random network. In the random network model of Beeman 
et. al. (|1984p , the structure consists of sp^ and sp^ sites in 
different proportions (see, also Robertson 1991). Mennella 
et al. ()1998p used continuous distributions of ellipsoids. 

Still another kind of disorder is obtained with pure 
PCG: this is orientation disorder. Papoular and Papoular 
()2009p showed that industry-made pellets of compressed 
powder of pure graphite (which can be considered as ar- 
tificially made disordered graphite) exhibit a feature at 
2200 1(4.52 /im"^), 1.3 /im"^ in width. This is one of 
the few largest widths observed in the sky by Fitzpatrick 
and Massa ([T986|) . 

It appears, therefore, that a vast gamut of carbonaceous 
materials display the UV bump, and that the latter's pro- 
file depends very much on the amount of structural and 
chemical disorder in the carrier. The problem, then, is to 
model this disorder as generally as possible, based on the 
best available solid state theory and in such a way as to al- 



low a straightforward derivation of the optical properties as 
a function of type and degree of disorder, at least over the 
whole range of the tt resonance, with the minimum possible 
ad hoc assumptions and adjustable parameters. This should 
also help understanding why most of the UV bump widths 
flock within a small interval and why they are limited from 
below. 

For this purpose, we start, here, with the carrier model 
consisting of stacks of graphitic layers in the form of ran- 
domly oriented small spherical grains, as studied by Draine 
and Lee ()1984p . These authors invoke dielectric functions 
derived from measurements on bulk graphite. In this paper, 
we derive these functions from an accurate tight-binding de- 
scription, based on the dramatic progress made in recent 
years in the theoretical and experimental study of graphene 
in single and multiple layers (SLG, MLG), the constitutive 
elements of graphite. Three of us (S.Y., R.R. and M.K.) 
have been engaged in this research for some time and have 
recently completed a theoretical model of the electronic and 
optical properties of SLG and MLG from the infrared to the 
UV, encompassing the tt resonance of sp^ bondings (Yuan et 
al. 2011a, 2011b). This model has been successfully applied 
to account for different experimental measurements, as the 
TT-plasmon resonance observed by electron energy loss spec- 
troscopy (EELS) (Eberlein et al. 2008 and Cass et al. 2008), 
the effect of disorder on the optical conductivity spectrum 
of doped graphane (Li et al. 2008), or the screening prop- 
erties and the polarization function obtained from inelastic 
X-ray scattering (Reed et al. 2010). The method allows for 
the consideration of different types of non-correlated and 
correlated disorder, such as random or Gaussian potentials, 
random or Gaussian nearest-neighbor hopping parameters, 
randomly distributed vacancies or their clusters, and ran- 
domly adsorbed hydrogen atoms or their clusters. We notice 
that, for the same amount of disorder, the electronic proper- 
ties of SLG and MLG are less affected by the Gaussian cor- 
related disorder, as the symmetry of the clean sample is less 
broken as compared to the non-correlated disorder (Yuan et 
al. 2011b). In fact, completely randomly distributed (non- 
correlated) impurities lead to the largest effect on the optical 
spectrum. Therefore, we consider here four different sources 
of non-correlated disorder, each of them controlled by a sin- 
gle parameter: 1) random on-site potential, corresponding 
to a local shift of chemical potential (inhomogeneous dis- 
tribution of electron and hole puddles), 2) random hopping 
parameter between two nearest neighbors, corresponding to 
random variations of inter-atomic distances and directions, 
3) randomly distributed vacancies, corresponding to atomic 
defects with missing carbon atoms (presence of sp^ bonds in 
the sp^ networks), and 4) randomly distributed hydrogen- 
like resonant impurities (including organic molecules such as 
CH3, C2H5, CH2OH, as well as H and OH groups). 

We consider mostly sp^-conjugated carbon in stacks 
of finite graphene sheets, arranged in the graphitic stack- 
ing schemes ABA (Bernal type), or ABC (rhombohedral) , 
which are sketched in Fig. [1] and perturbed by the differ- 
ent kinds of disorder introduced above. While our results 
clarify a number of issues concerning graphene and graphite 
properly, here we are mainly interested in the impact of im- 
purities and disorder on the position and width of the tt 
surface plasmon (UV bump). As is well known (Katsnelson 
2012), the band structure of graphene presents saddle points 
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Figure 1. Atomic structure of the two different stacking sequences of the graphene layers considered in this work: ABA- (left) and 
ABC-stacked (right) MLG. The nearest neighbor intra-layer t and the inter-layer 71 and 73 hopping amplitudes are schematically shown 
in the figure. 



(Van Hove singularities), with the corresponding logarith- 
mic divergencies in the density of states (DOS), at energies 
E — ^t (in the nearest neighbor approximation), where t 
is the nearest neighbor hopping parameter. The 7r-plasmon 
mode is associated to electron-hole transitions between the 
Van Hove singularities of the valence (tt) and the conduction 
(tt*) bands. As the disorder parameter increases, it is found 
that the slope of the blue wing of the tt resonance decreases 
in such a way as to increase the width of the plasmon fea- 
ture, leaving its peak position nearly unchanged, as required 
by astronomical observations. This trend is common to all 
types of disorder considered, but is more or less pronounced. 
Our results show that random on-site potential disorder is 
the most effective, but, for all types, at sufficiently high de- 
grees of disorder, the tt resonance is completely damped, 
indicating that the range of our computations covers the 
most extreme situations, from perfect order to very high 
disorder observed on laboratory materials. The UV bump 
survives considerable disorder, at the price of broadening 
and frequency shift. 



This paper is organized as follows. Section 2 lays down 
the basic principles, assumptions and approximations under- 
lying our calculations, and discusses various types of disor- 
der (imperfections) considered here. Section 3 displays some 
of our results in the form of the electrical conductance of 
graphene as a function of normalized wave number, for dif- 
ferent types and degrees of disorder. Section 4 describes the 
procedure which delivers the corresponding dielectric func- 
tion and, hence, the extinction cross-section. Section 5 dis- 
plays a number of computed UV bump profiles and high- 
lights their changes as a function of disorder. In Sec. 6, our 
computed results are used to derive the extinction bumps 
of free- flying graphitic bricks, and these are compared with 
those of a few disordered laboratory model carrier candi- 
dates. Finally, in Sec. 7, all these cases are confronted with 
available astronomical observations. A discussion of these re- 
sults and their possible improvements and applications con- 
cludes this work . 



2 THEORETICAL PROCEDURE 

We consider the effect of different sources of disorder on SLG 
and MLG by direct numerical simulations of electrons on a 
honeycomb lattice in the framework of the tight-binding ap- 
proximation. This allows us, by means of the time-evolution 
method, to obtain the DOS of large samples containing mil- 
lions of atoms. The time-evolution method is based on the 
numerical solution of the time-dependent Schrodinger equa- 
tion (TDSE) with additional averaging over random super- 
position of basis states (Hams et al. 2000, Yuan et al. 2010). 
The electronic single particle dispersion in SLG is properly 
described by the hopping of electrons between nearest neigh- 
bor carbon sites arranged in a honeycomb lattice with two 
atoms per unit cell. For MLG, additional inter- layer hop- 
ping processes must be included in the kinetic Hamiltonian. 
A sketch of a MLG with the different hopping parameters 
used in the calculation is shown in Fig. [T] 

The tight-binding Hamiltonian of a MLG is given by 

H^ Y. Hi+ Y. H[ + H,^^, (1) 

1^1 1=1 

where Hi is the Hamiltonian of the /'th layer of graphene, 

Hi^- ^(ti,i^-a[.6ij + h.c) + ^^;i,iq^Q,i, (2) 

where a\ ^ {bi,i) creates (annihilates) an electron on sublat- 
tice A (B) of the /'th layer, and ti^ij is the nearest neighbor 
hopping parameter if i, j are nearest neighbors and zero oth- 
erwise. It is known that the nearest-neighbor approximation 
is quite accurate to describe electronic structure of singla- 
layer graphene (Castro-Neto et al. 2009, Katsnelson 2012). 
The second term of Hi accounts for the effect of an on-site 
potential vi^i, where ni^i — c\ ^ci^i is the occupation number 
operator and c = a, 6. In the second term of the Hamilto- 
nian Eq. ([1]) , H'l describes the hopping of electrons between 
layers / and / + 1. We consider the two more stable forms of 
stacking sequence in bulk graphite, namely ABA (Bernal) 
and ABC (rhombohedral) stacking, as shown in Fig. [1] For 
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a MLG with an ABA stacking, H'l is given by 

3 3,3' 

(3) 
where the inter- layer hopping terms 71 and 73 are shown 
in Fig. [U Thus, all the even layers (/ + 1) are rotated with 
respect to the odd layers (/) by +60°. The difference between 
ABA and ABC stacking is that, the third layer (s) is rotated 
with respect to the second layer by —60° (then it will be 
exactly under the first layer) in ABA stacking, but by +60° 
in ABC stacking (Yuan et al. 2011a). The spin degree of 
freedom contributes only through a degeneracy factor and 
is omitted for simplicity in Eq. ([T]). 

The term Himp describes the hydrogen-like resonant im- 
purities: 



De Raedt 2000, Yuan et al. 2010) 



H^mp = ed^ d\d^ + ^ ^ {d\c^ + h.c) , 



(4) 



where Sd is the on-site potential on the impurity and V 
is the hopping between carbon and absorbed atom. The 
band parameters V ^ 2t and ea ~ — 1/16 are obtained 
from the ah initio density functional theory (DFT) calcu- 
lations (Wehling et al. 2010). Following Refs. (Wehhng et 
al. 2010, Yuan et al. 2010, 2011b), we call these impurities 
as adsorbate hydrogen atoms but actually, the parameters 
for organic groups are almost the same (Wehling et al. 2010). 
The band parameters for OH-group and fluorine are similar 
to those for hydrogen except that the on-site potential on 
impurities is td ~ —t (Yuan et al. 2012). A vacancy is an 
atomic defect with one carbon atom missing, and can also 
be regarded as an atom (lattice point) with on-site energy 
Vi ^ 00 or with its hopping parameters to other sites being 
zero. 

The numerical calculations of the optical conductivity 
and DOS are performed based on the numerical solution of 
the TDSE for the non-interacting electrons. The frequency- 
dependent optical conductivity follows (Ishihara 1971, Yuan 
et al. 2010) 
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(5) 



(we take h = 1) where /S = l/ksT is the inverse tempera- 
ture, Q is the sample area, / {H) = 1/ L^^^-^) + ij is the 

Fermi-Dirac distribution operator, J a (t) — e^^^J^e"*^* is 
the time-dependent current operator in the a {— x 01 y) 
direction, and \Lp) is a random superposition of all the ba- 
sis states in the real space, i.e., (Hams and De Raedt 2000, 
Yuan et al. 2010) 



\^) = Y.''^C\^., 



(6) 



where ai are random complex numbers normalized as 
5^^ \ai\ =1, and |0) is the vacuum state. The time evolution 
operator e~*^* and the Fermi-Dirac distribution operator 
/ (if) are obtained by the standard Chebyshev polynomial 
representation (Yuan et al. 2010). 

The density of states is calculated by the Fourier trans- 
form of the time-dependent correlation functions (Hams and 



1 C" 



' {if\e '"^\<fi)dt, 



(7) 



with the same initial state \(p) defined in Eq. JB}. For a more 
detailed description and discussion of our numerical method 
we refer to Ref. (Yuan et al. 2010). In this paper, we fix the 
temperature to T = 300K. In our numerical calculations, 
we use periodic boundary conditions in the plane (XY) of 
graphene layers, and open boundary conditions in the stack- 
ing direction (Z). The size of SLG in our calculations is 
8192 X 8192 or 4096 x 4096. For MLG, we keep the total 
number of atoms in the same order as for SLG. 

Each kind of disorder considered here is controlled by a 
single parameter: 1) random on-site potential is controlled 
by Vr'. the potential Vi on each site can randomly change 
within the range [—Vr^Vr]; 2) random hopping tr'. the hop- 
ping between two nearest neighbors can randomly vary 
within the range [t — tr,t -\- tr]; 3) concentration of random 
distributed vacancies Ux = Nx/N, which is the probability 
for an atom to be replaced by a vacancy in the lattice, where 
Nx is the total number of vacancies and N is the total num- 
ber of sites in the lattice; 4) similarly, the concentration of 
random distributed hydrogen ad atoms, is controlled by the 
parameter rii = Ni/N, which is the probability for an atom 
to absorb a hydrogen atom. 

It should be stressed that our calculations automati- 
cally account for the variation of the matrix element of the 
dipole moment, which is often disregarded for simplicity (see 
Bassani F. and Pastori-Parravicini G. 1967; Mennella et al. 
1998), but is essential for our present purposes. 



3 COMPUTATIONAL RESULTS: 
CONDUCTANCE 

The most direct test of the theory is the measurement of the 
absorbance of light perpendicularly through the graphene 
layer, which is proportional to the conductance of the latter. 
For this aim, in this section we present numerical results for 
the optical conductivity, obtained from Eq. (|5]), for SLG and 
MLG with different kinds of disorder. Our results are shown 
in Fig. [2 

We first consider the effect of random local change of 
on-site potentials and random renormalization of the hop- 
ping, which correspond to the diagonal and off-diagonal dis- 
orders in the single-layer Hamiltonian Eq. [2l respectively. 
The former acts as a local shift of the chemical potential of 
the Dirac fermions, and the latter is associated to changes 
of distance or angles between the pz orbit als. As we have ex- 
plained above, we introduce non-correlated disorders in the 
on-site potentials by considering that the on-site potential 
Vi is random and uniformly distributed (independently of 
each site i) between the values —Vr and -\-Vr. Similarly, non- 
correlated disorder in the nearest-neighbor hopping is intro- 
duced by letting tij be random and uniformly distributed 
(independently of the pair of neighboring sites {i,j)) be- 
tween t — tr and t-\-tr. Both kinds of disorder have a rather 
similar effect on the optical conductivity, as it can be seen 
from Figs. [2ja) and (b). The spectrum is smeared at the 
resonance energy cj = 2t, and the smeared region expands 
around the singularities neighboring areas as the strength 



Effects of disorder on grain spectra 5 







^ 
^ 




Figure 2. Optical conductivity cr(u) of SLG with different kinds of disorder, in units of ctq = 7re^/2/i. (a) Effects of random on-site 
potentials on (j(uj). The on-site potential is let to randomly vary in the range [—Vr,Vr], for Vr = 0, 0.2, 0.3, 0.4, 0.5, 0.7, 1, 2.5 (in units of 
t). (b) Effects of random hopping constants on a(uj). The hopping integral is let to vary randomly between [t — tr,t-\- tr], for tr = 0, 0.3, 
0.35, 0.4, 0.45, 0.5 (in units of t). (c) Effect of random vacancies (missing atoms) in the lattice, with the concentration per carbon atom 
Ux taking the values 0, 0.01, 0.05, 0.1. (d) Effects of random hydrogen-like impurities (adatoms or admolecules) with the concentration 
per carbon atom rii being 0, 0.01, 0.05, 0.1, 0.2, 0.5, 1.0. For all the plots, the value of the corresponding disorder parameter Vr, tr, fix 
and rii increases in order of decreasing sharpness of the peak at lj = 2t. 



of the disorder increases. We remind here that the origin of 
this resonance is associated to particle-hole excitations be- 
tween states of the valence band with energy E ^ —t and 
states of the conduction band with energy E ^ t, which 
contribute to cr(Lu) with a strong spectral weight due to the 
enhanced density of states at the Van Hove singularities of 
the TV bands, in agreement with recent experimental results 
(Mak et al. 2008). 

We next consider the influence of vacancies and hydro- 
gen impurities on the spectrum. Apart from the creation 
of mid-gap states in the DOS, which is not relevant for the 
present discussion, the presence of vacancies and/or adatoms 
leads also to a smearing of the Van Hove singularities in the 
DOS. This is in fact the behavior found in Fig.[2jc) and (d), 
where we show the optical conductivity of SLG with differ- 
ent concentrations of impurities. In particular, the presence 
of low or moderate concentrations of hydrogen impurities, 



which are introduced by the formation of a chemical bond 
between a carbon atom from the graphene sheet and a car- 
bon, oxygen, or hydrogen atom from an adsorbed organic 
molecule (CH3, C2H5, CH2OH, as well as H and OH groups), 
have a quite similar effect on the electronic structure and 
transport properties of graphene as compared to the effect 
of vacancies. In our work, the adsorbates are described by 
the Hamiltonian Himp, Eq. Q, where the band parameters 
V ^ 2t and Sd ~ — 1/16 are obtained from ab-initio density 
functional theory (DFT) calculations (Wehling et al. 2010). 
On the other hand, large concentrations of hydrogen impu- 
rities change dramatically the optical properties of the sys- 
tem. In fact, for the maximum concentration rii = 1, which 
correspond to an H atom bonded to each C atom of the lat- 
tice, we have no longer graphene but a completely different 
material, called graphane which contrary to graphene, has 
a band gap of A ^ 2t. Therefore, the optical conductivity 
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Figure 3. Effect on (t{ijo) of the number of layers and stacking order of clean multi-layer graphene. Plots (a) and (b) are for the ABA 
stacking, whereas plots (c) and (d) are for ABC. Plot (b) and (d) are zooms of the spectra (a) and (c), respectively, around the resonance 
peak at cj ~ 2t. The main effect, for present purposes, is a local deformation of the peaks [zooms (b) and (d)]. The blue wing, on which 
the Frohlich resonance stands, is hardly affected. The maximum separation between sub-peaks of the resonance top is determined by the 
coupling energy between adjacent layers, 71 = 2t/15. 



at low frequencies is zero from cj = to cj = 2t, and from 
this energy it grows smoothly, presenting a series of peaks 
associated to newly created Van Hove singularities in the 
spectrum, as it can be seen by looking at the purple line of 
Fig. [2jd). However, the more realistic case of low or moder- 
ate concentrations of hydrogen, leads to a smearing of the 
Van Hove singularity with a transfer of spectral weight to 
new peaks in the spectrum, particularly to a new peak at 
cj ~ t which is associated to transitions between the im- 
purity zero energy band and the Van Hove singularities at 

As an infinite sheet of graphene proper cannot survive 
in vacuum (it folds over itself), we have explored the elec- 
trical effects of stacking an increasing number of sheets. In 
equilibrium, two stacking schemes must be considered, ac- 
cording to the arrangement of the benzenic rings facing each 
other: ABA and ABC (see Fig. [Tj). The inter- layer distance, 
experimentally measured, has a value oi d ^ 3.36 A. The 
inter-layer hoppings considered here are 71 = 2t/15 and 



73 = t/10^ the same values as in graphite, and are also 
sketched in Fig. [T] The stacking order, which leads to differ- 
ent low energy properties of ABA and ABC MLGs, has very 
limited effect on the high energy spectrum and, in particu- 
lar, on the blue wing of the tt resonance, where the Frohlich 
bump stands. This can be seen in Fig. [3ja) and (c), which 
show the full oiuo) spectrum of MLG with an increasing 
number of layers, up to 100 layers, with ABA- and ABC- 
stacking order, respectively. The UV bump, therefore, will 
not be affected. By contrast, the effect on the peak is re- 
markable, and is a classical outcome of coupling identical 
oscillators: creating an increasing number of subpeaks on 
both sides of the single-oscillator peak. Note, however, that 
their excursion from the latter does not exceed an amount 
determined by the coupling energy assumed between lay- 
ers (see Griineis et al. 2008), ^ 3% in the present case. As 
a consequence, the peak becomes less sharp. This must be 
a major reason why the tt and a resonances of graphite are 
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much less peaked than those of graphene (see also Kobayashi 
and Uemura 1968; Johnson and Dresselhaus 1972). 

Summarizing, clean graphene has the steepest slope of 
the TT-feature. This slope is shown below to govern the width 
of the Frohlich resonance; so, we will concentrate on those 
types of disorder which have the strongest effect on the slope 
of the blue wing. The most striking aspect of these curves, 
for our present purposes, is the continuous variation of the 
slope of the high energy wing of the tv resonance, as the dis- 
order parameter increases. All types of disorder exhibit this 
behavior, albeit to different extents. Other, so-called corre- 
lated disorders were also studied, but they do not affect the 
blue wing notably. The reason is that for correlated disorder, 
for which the distribution of the disorder follows particular 
topological structures, such as Gaussian potentials or Gaus- 
sian hopping parameters, or resonant clusters with groups of 
vacancies or hydrogen adatoms, the electronic DOS of clean 
graphene is highly preserved. Therefore, a SLG or a MLG 
sample with a given amount of impurities, will present an 
optical spectrum closer to that of clean graphene (in partic- 
ular presenting a prominent peak at the vr-resonance) if the 
impurities are correlated, rather than if they are spread in 
a non-correlated manner in the sample. Another remarkable 
effect of disorder is the progressive damping of the tt reso- 
nance itself, which is a very well experimentally documented 
phenomenon exhibited by families of a-C:H materials. We 
emphasize that no ad hoc broadening mechanism, such as 
an arbitrary life time, is invoked in the present treatment. 
Each single parameter controlling any given type of disor- 
der can be estimated from the structural and environmental 
characteristics of the corresponding sample. 



4 THE DIELECTRIC FUNCTION 

The study of the Frohlich resonance requires the knowledge 
of the dielectric function. The imaginary part of the relative 
dielectric function, e = ei+ie2 is related to the conductance 
computed above by 



a 



(8) 



where a and u are in units of s~^ in the SI system. While 
the conductance a is completely defined in terms of ao = 
7re^/2/i, which depends only on fundamental constants, the 
corresponding energy variable is given only to a factor t, 
which represents the coupling energy between neighboring 
atoms in perfect, clean graphene. The latter is not known 
with great accuracy, and the value adopted for it differs ac- 
cording to authors (see, for instance, Griineis et al. 2008, 
Kobayashi J. and Uemura Y. 1968). Here, we take 2t = 4.21 
eV, so that the tt peak of the dielectric function falls near 
the measured position for graphite (Taft and Philipp 1965; 
Greenaway et al. 1969; Tosatti and Bassani 1970). 

The higher energy a resonance of graphite peaks near 
14 eV but its red wing slightly overlaps the blue wing of 
the TT resonance, in between them, and this bears upon the 
position of the Frohlich resonance. However, the theoreti- 
cal procedure outlined in Sec. 2 has not yet been extended 
to the a resonance (see Sec. 7.7). We were therefore led to 
use, instead, the values of 62 deduced by Taft and Philipp 
()1965p from reflectance measurements between 9 and 32 eV, 
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Figure 4. Effects of random on-site potentials on the dielectric 
function. The disorder parameter, Vr (in units oft), is 0, 0.2, 0.3, 
0.4, 0.5, 0.7, 1, 1.5 2.5, in order of decreasing sharpness of the 
resonance peak. 



and joined them continuously to each of the tt curves we 
computed as explained above. By extending the overall fre- 
quency range, this also helps improving the accuracy of the 
Kramers-Kronig transformation which was applied to 62 in 
order to deduce ei. 

This anzatz is justified by the fact that the intensity of 
the a resonance of carbonaceous materials is largely decou- 
pled from that of the tt resonance, as the two singularities 
are separated by more than 6 eV; a fact which is likely to 
be related with the lack of correlation, in the interstellar 
extinction curves, between the UV bump and the far UV 
curvature. Moreover, disorder is not likely to perturb exces- 
sively the a resonance since, as recalled in the Introduction, 
the latter is associated with the ppa type of bonding which 
survives the transition to sp^ hybridization that dominates 
structurally disordered materials. This statement is borne 
out by the comparison of the VUV spectra of graphite with 
several species of coals of different evolutionary states, which 
correspond to different degrees of disorder Papoular et al. 
1995). 

Figures 4 to 6 exhibit some of the results for 3 types of 
disorder, excluding the vacancy disorder, which was shown 
above not to affect significantly the slope of the blue wing of 
the conductance resonance. They mirror the curves of Fig. 
[2] (a), (b) and (d). The general trends of interest to us in the 
three graphs are the same, even though their magnitudes 
may differ: the peak for pure and perfect graphene is a loga- 
rithmic singularity; as disorder increases, it is progressively 
damped, and ultimately drowns into a continuum; concur- 
rently, the blue slope becomes less and less steep. The former 
behavior may be ascribed to the progressive suppression of 
sp^ bonds, and accounts for the observed difference between 
graphite and various species of a-C:H (see, Robertson 1986, 
Mennella et al. 1998) or coals (see Papoular et al. 1995). The 
change in slope is shown below to account for the width vari- 
ations of the Frohlich resonance feature. Note that this effect 
is all but absent in the case of vacancy disorder. 
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Figure 5. Effects of random hopping constants on the dielec- 
tric function. The disorder parameter, tr (in units of t), is 0.2, 
0.3, 0.35, 0.4, 0.45, 0.5, in order of decreasing sharpness of the 
resonance peak. 



Figure 7. Effects of random on-site potentials on the Frohlich 
resonance of the tt resonance. The disorder parameter, Vr (in units 
oft), is 0, 0.4, 0.7, 1, 1.5, 2.5, in order of decreasing sharpness of 
the resonance peak. 
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Figure 6. Effects of random atomic impurities on the dielectric 
function. The disorder parameter, n^, is 0.01, 0.05, 0.1, 0.2, in 
order of decreasing sharpness of the resonance peak. 



5 APPLICATION TO THE 2175-A FEATURE 

The weak scattering apparently included in the astronom- 
ical UV bump extinction, elicits the assumption of small 
particles for the band carrier, and, hence, the application of 
the Rayleigh approximation. Consequently, surface effects 
cannot be neglected, which entails the use of the Frohlich 
expression for the extinction efficiency of a small sphere of 
radius a: 

rM 247r 62 .^x 

^/«=— (e,+L--l)^ + M^ (9) 

where A is the wavelength of the radiation, ei,2 are, respec- 
tively, the real and imaginary parts of the dielectric function 
of the assumed homogeneous material; L is the shape param- 
eter of the supposedly ellipsoidal grain, 1/3 in the case of 
a sphere (see Bohren and Huffman 1983); a is the radius 
of the sphere. In this approximation, there is no need to 
distinguish between absorption and extinction efficiencies, 
both being designated by Q. The real part, ei, of the dielec- 



tric function is obtained by applying the Kramers-Kronig 
transformation to 62 from to 32 eV. For this operation, 
Papoular (|2009p built a transform algorithm based on the 
Maximum Entropy method. Figures 7 to 9 show the results 
of computations for a few representative cases. The "clean" 
case (all disorder parameters null) is represented in Fig. [71 
only (black line); its peak is slightly higher than that for 
m = 0.01 in Fig. [9] 

As expected, as disorder increases, the width of 
the feature increases and the height decreases, both 
monotonously; but the peak position hardly changes from 
~ 4.6 //m~^ (because of the Kramers-Kronig relation be- 
tween ei and 62). The minimum width and maximum height 
are carried by the clean case. For this case, the peak value of 
Q/a is ^ 0.02 A~^ (upper limit of Q/a in our study) and the 
width, ^ 0.85 /im~^ (lower limit in the present work). These 
limits are set by the logarithmic slope of the resonance (see 
van Hove 1953) . 

Strictly speaking, these results apply to a single 
graphene layer. However, they derive solely from the behav- 
ior of the blue wing of the conductance resonance, and this 
was shown to be insensitive to the number of stacked layers 
[Fig. [3] (a) and (c)], whatever the stacking scheme. Figure 
10 illustrates this point. We therefore assume, below, that 
they apply equally to MLG (multi- layer graphene). 



6 APPLICATION TO ASTRONOMY 

Consider first the free- flying brick model (Fig. [6]), where each 
dust particle is a spheroidal assembly of nano-sized portions 
of graphene, stacked according to the ABA or ABC scheme. 
Their orientation being random, the apparent extinction 
along a given line of sight will be the average of extinc- 
tion for E ^- c (electric field parallel to the layer planes) and 
for E\\c (electric field perpendicular to these planes), with 
relative weights 2/3 and 1/3, respectively. We have, there- 
fore, applied the same theoretical procedure to compute the 
dielectric function in the parallel case. As an example, Fig. 
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Figure 8. Effects of random hopping constants on the Frohhch 
resonance of the tt resonance. The disorder parameter, tr (in units 
of t), is 0.2, 0.3, 0.35, 0.4, 0.45, 0.5, in order of decreasing sharp- 
ness of the resonance peak. 
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Figure 9. Effects of random atomic impurities on the Frohhch 
resonance of the tt resonance. The disorder parameter, Ui, is 0.01, 
0.05, 0.1, 0.2, in order of decreasing sharpness of the resonance 
peak. 



Figure 10. The Frohhch resonance of the tt resonance in the 
clean case for a single layer (red), a bi-layer (green) and an ABA 
stack of 100 layers (blue). 





Figure 11. Models of dust particles. A: free- flying bricks ; each 
dust particle is made of one ordered stack of graphene layers. 
B: agglomerated bricks; each dust particle includes several such 
"bricks" . 



2009) . It is made of powder of nano-sized highly graphitized 
and pure carbon grains, compressed under high pressure. 
The measured reflectance of PCG was modeled by 2 pairs of 
Lorentzian oscillators, for the parallel and perpendicular tt 
and a dielectric functions. The isotropic dielectric function 
is a mix of the parallel and perpendicular dielectric func- 
tions, and was computed using the Bruggeman mixing for- 
mula. The extinction feature was then found to peak at 4.52 
/im"\ with FWHM 1.3 /xm"^ ("PCG agglom" in Tab. 1). 



[12] plots the real and imaginary parts of this function for the 
clean case. 

Correspondingly, Fig. [13] plots the parallel and perpen- 
dicular extinction efficiencies for the clean case, together 
with their weighted average. Obviously the position and rel- 
ative weakness of the "parallel bump" are such that the 
average curve differs from the "perpendicular bump" only 
in that its height is reduced by a factor 2/3. For purposes of 
quantitative comparison, the average Q/a was fitted with a 
"Drude" curve, as initiated by Fitzpatrick and Massa (|1986p , 
also shown in the figure. It peaks at 4.65 /xm"^ and its width 
is 7 '^ 0.85 /im~^. Only for highly disordered structures, 
when the "perpendicular" bump is weakened beyond the 
range considered in Fig. 7-9, does the effect of the "parallel" 
bump on the average become notable. 

In another dust model (Fig. IIB), the bricks are ag- 
glomerated into a single particle: industrial poly crystalline 
graphite (PCG) is a case in point (Papoular and Papoular 
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Figure 12. Dielectric function for the clean case, and for E\\c. 
Red line: imaginary part 62; green line: real part ei. Color on line. 
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Figure 13. Blue line: average extinction bump for randomly ori- 
ented spheroidal stacks of clean graphene layers. Black points: 
Drude fit. Red upper line: extinction bump for E in plane. Green 
lower line: extinction bump for E normal to the plane. Color on 
line. 



Table 1. Experimental model results. 



Name 


z^o 


7 


AR^ 


4.46 


1.48 


BCC V ^ 


4.7 


1.6 


Free bricks^ 


4.56 


1.15 


PCG agglom^ 


4.52 


1.3 



The units of uq and 7 above are iim ^ . 
a Schnaiter et al. (1998]) 
b Mennella et al. (1998)) 
c,d Papoular and Papoular (|2009)) 



Table 1 collects examples of relevant experimental results, 
also symbolized in Fig. 14. 

It is instructive to envision the behavior of this material 
if each of its constitutive graphitic bricks were set free in the 
atmosphere, as in model 10 A. The parallel and perpendic- 
ular extinctions in that case were calculated using the same 
dielectric functions as deduced above from the reflectance 
of PCG. Their weighted average (1/3, 2/3) was found to 
exhibit a bump at 4.56 /im~\ with FWHM 1.15 /im~\ Ob- 
viously, agglomeration of the bricks entails some broadening 
as do the various types of disorder considered above. There 
is also some redshift. Unlike with the free- flying bricks, the 
effective dielectric function, here, is very sensitive to the par- 
allel dielectric functions. 

A similar behavior is displayed by nano-sized HAC par- 
ticles agglomerated in pure argon atmosphere (Schnaiter et 
al. 1998): the feature width increases with the degree of 
clustering (their Fig. 3; the narrowest of these features is 
included as"AR" in Tab. 1 and symbolized in Fig. 14). Ag- 
glomerates obtained with hydrogen included in the atmo- 
sphere show less clear trends. 



7 DISCUSSION 

7.1 Comparison with observations 

All cases computed above are represented in Fig.[T4]as points 
with coordinates (7,zyo). For comparison, the points plotted 
by Fitzpatrick and Massa (|2007p in their Fig. 16 are also 
represented (schematically). While there is a rough agree- 
ment in the width distribution (abscissae), the ordinates, 
z^o, of the computed points are systematically shifted up- 
wards by ^ 0.05 A, or about 1 %. z/q depends slightly on 
the fundamental hopping parameter t (Sec. 2), and mostly 
on the electron plasma frequency of the system, which, in 
turn, depends on the shape of the computed electronic bands 
through the effective number of electrons, Ueff (see Taft and 
Philipp 1965). In the tight binding scheme adopted here, this 
is found to be 0.56 electron/atom, a value which is compa- 
rable to that determined experimentally for graphite, by ap- 
plying the e2-Sum Rule (see Altarelli et al. 1972). According 
to the latter, Ueff is linearly related to 62- 

Calculation shows that uq can be brought down to the 
right value by taking Ueff = 0.53, which amounts to re- 
ducing 62 uniformly by 6 %. Although this is well within 
experimental errors of measurement of €2 , we chose to keep 
Ueff = 0.56, so as to let one assess the validity of the tight 
binding assumption. As for the minimum width determined 
here, 0.85 /xm"^, we note that it is slightly higher than the 
value ^ 0.8 /im~"^ suggested by astronomical measurements. 
This is perhaps one weakness of our model, or of the tight- 
binding approximation used in the derivation. 

On the other hand, it should be stressed that, apart 
from the selected disorder type and parameter, which define 
the bump carrier structure, the only adjustable parameter 
in the algorithm is the nearest-neighbor hopping parameter, 
t (half the transition energy at resonance) . Even this is re- 
stricted by experiment within tight limits. This ensures that 
the relation between disorder and bump properties remains 
clear: disorder damps the resonance, which flattens its high- 
energy wing; hence, the UV bump is weakened and broad- 
ened. The reason why the peak position does not change 
much lies in the Kramers-Kronig relation between a and 

€2- 

Figure [Ml shows that, for moderate degrees of disorder, 
the representative points flock near the clean case (minimum 
bump width) , in about the same width range as the observed 
widths. This sets upper limits to the extent of disorder sug- 
gested by the present model for bump carriers. Estimates of 
these limits for the parameters of the three types of disorder 
considered here are, respectively: Vr ^ 0.7t, tr ^ 0.3t and 
n^ - 7%. 

For higher values of Vr (not shown in Fig. I14p it is 
found that the peak height and wavenumber decrease, and 
the width increases, with disorder. This is in line with the 
measurements on coals in order of decreasing graphitization 
(see Papoular et al. 1993). Figures 4 to 6 for 62 also suggest 
that highly disordered grains provide a UV continuum, as 
expected for strong resonance damping, also in accordance 
with the measurements on coals. Such grains may contribute 
to the ISEC continuum. 

Mennella et al. ()1998p obtained HAC samples from elec- 
trical discharges and measured their extinction efficiency, 
from which they deduced the dielectric function and the 
Frohlich resonance. Even the vr resonance hardly emerges 
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Figure 14. Peak position vs FWHM of the UV bump. Observa- 
tions (adapted from Fig. 16 of Fitzpatrick and Massa 2007): the 
red eUipse encloses the vast majority of observed representative 
points; red crosses: a few out hers. Our computations: Open cir- 
cle: randomly oriented clean and ordered graphene bricks. Filled 
squares and solid line: random on-site potentials (vr = 0, 0.2, 
0.4, 0.5, 0.7, 1, 1.5, 2.5, 5); filled circles and dotted fine: random 
hopping parameters (tr = 0.1, 0.2, 0.3, 0.35, 0.4, 0.45, 0.5); filled 
triangles and dashed line: random foreign atoms (e.g. hydrogen; 
ni=l, 5, 10, 20 %). Laboratory models: open triangles (see Table 
1). Color on line. 



in nascent samples. After thermal annealing, it is distinctly 
visible, and the UV bump shows up. Under further treat- 
ment by UV irradiation, the latter becomes narrower and 
stronger; this is included in Tab. 1 as BCC V and symbol- 
ized in Fig. [21 In order to obtain still narrower features 
the authors extrapolated the dielectric function as a func- 
tion of radiation dose, and the correspondingly computed 
widths were then compatible with observed bump widths, 
with only a slight decrease in peak wavenumber. This trend 
is analogous to our findings. It is attributed to an increase in 
aromatic cluster number and average size, and a reduction 
of size dispersion, which we interpret as reduced disorder. 

Several other causes of bump broadening have been con- 
sidered in the past. Sorrell (|199Qp , for instance, argued that 
electron-phonon scattering ("viscous flow") and collision of 
electrons with adsorbed H atoms limited the bump width 
from below. The latter phenomenon is a close kin to our im- 
purity disorder, while the former is not considered here for 
we are only interested in hopping electrons, not free (con- 
duction) ones. Another cause may be the smallness of the 
grain size. 

The argument developed above hinges upon the 
graphene properties. For present purposes, this was also 
shown above to apply to stacks of single- layered graphene 
(MLG). However, graphene theory assumes an infinite plane, 
while, for our extinction calculations to apply, the bump car- 
rier size is limited to ^ 100 A by the Rayleigh condition 
for 2175 A . In fact, it can be shown that the minimum size 
of a graphene layer for the properties of an ideal layer to be 
recovered is ~ 15 A, the size of a big PAH (see Sec. 7.3). 

While other causes may coexist in most astronomical 
instances, our model excludes their contribution in cases 
corresponding to the minimum observed bump widths, as 
the clean case already has this minimum width. 



7.2 Carbon budget 

We seek now to estimate the fraction of the available carbon 
that must be locked in this model bump carrier for it to ac- 
count for astronomical observations. Fitzpatrick and Massa 
([2D0 7) have cast their interstellar extinction measurements 
in the form of a normalized color excess (in units of magni- 
tude) 
k{X - V^) = E{B^V) ^ <^i +C2X + C3i^(x;7,Xo) + C4F(x), 

2 

where x = A~^, Ao is the bump peak wavelength, the c's are 
constant factors for each line of sight and F(x) represents the 
far UV curvature of the extinction curve. Admittedly, this 
mathematical representation implies no particular physical 
model of the "bump"; however it is of great help in com- 
paring quantitatively our model results with astronomical 
observations. Since our model carries a negligible amount 
of underlying continuum (see Fig. 7 to 9), it can only con- 
tribute to the term C3D. The peak value of this term, at Ao, 
is 

AA:(Ao -V) = C3/7' 

and corresponds to the peak value of Q/a. Now, the Lilley 
empirical relation gives 
EiB-V) = ^ 

on average through the ISM, with N the line density of H 
atoms in cm~^. Hence, 



A^(Ao -V) = Ak{Xo - V) E{B - V) : 



C3 Nh 

72 51021 



(10) 



Also note that, by definition, 

^(A -V) = 1.09 (r(A) - T(y)), 

where r is the optical depth. So, AE{X — V) = 1.09 Ar, 

where At is the bump in the optical depth. 

On the other hand, if Ngr is the line density of bump 
carriers along the line of sight, and if Rayleigh's approxima- 
tion applies, then the bump in the optical depth is 
At = Tia^ NgrQ, 

assuming the grains are spherical with radius a, and there 
is no line saturation. Also, if p is the density of the grains, 
the line density of carbon atoms locked in the grains is 



Nc = Ng 



Ana^p 



3x 12 X 1.710-24 



(11) 



(12) 



Combining eq. 10 and 11, one finally obtains 

^C ^ C3P 

Nh l^{Q/ay 

One can see, from the compilation by Fitzpatrick and 
Massa ()2007p . that the average values of cs and 7 are, re- 
spectively, about 3 and 1 pm~^ (the corresponding bump 
area is 5/im~-^). Take 2 gm/cm^ for the density of graphite 
and, accounting for some disorder in our model grains, as- 
sume the average peak value of Q/a is 1.5 10^ cm~^. Then, 
on average, the fraction of cosmic carbon required to be in 
the model grains turns out to be ^^ 80 ppm, or about 1/4 
of the available carbon (for comparison, see Snow and Witt 
1995). Of course the largest fraction must be reserved for 
molecules and for the less ordered and pure carbon grains. 
As stated earlier, the latter retain their a resonance well be- 
yond the complete damping of their tt resonance. They may, 
therefore, well account for the observed rise of interstellar 
extinction towards shorter wavelengths. 
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Figure 15. Binding energy per atom (Kcal/mol) as a function of 
the number of benzenic carbon rings in planar, compact clusters; 
1 Kcal/mol = 0.0434 eV. 



7.3 Minimum graphene layer size 

The present model invokes stacks of small graphene sheets. 
We must, therefore, inquire whether the properties of 
graphene still apply to sizes smaller than the Rayleigh limit. 
These properties are determined by the interactions between 
an atom and its neighbours. Obviously, these interactions 
fade away as the distance to the central atom increases. In 
order to find the "cut-ofF' distance, we consider compact 
clusters of benzene rings of increasing size. With the help of 
a commercial Molecular Modeling package [Hypercube, Hy- 
perchem 7 (1996')], their structure was optimized by seeking 
the minimum total binding energy, using the Quantum Me- 
chanical semi-empirical AMl/UHF method. Figure 14 shows 
the variation of the binding energy per atom, e^, as a func- 
tion of the number of rings, Nr. 

Beyond about 40 rings, or 84 atoms, the binding en- 
ergy reaches, for all practical purposes, an asymptotic value 
of -182 Kcal/mol or 8 eV per atom, which is indeed about 
the graphite sublimation energy. The corresponding diame- 
ter of the structure is about 15 A, the size of a large PAH 
or Piatt particle (see Donn 1968). A similar conclusion was 
reached by Robertson (|1986p on the basis of a simpler Hiickel 
calculation (see his figures 22 and 28) . Thus the minimum 
extension required for a graphene sheet is small enough that 
many such sheets can be accommodated in a Rayleigh-sized 
grain. 



7.4 Effects of grain shape 

The Frohlich formula for Q/a includes the ellipsoidal charac- 
teristic parameter, L, which was set above at 1/3, assuming 
a spherical shape for both the grain and its subgrain com- 
ponents. As noted by Bohren and Huffman ([1983]), shape 
effects are no weaker in the Rayleigh approximation than 
for large particles, and the band shape of discs or needles 
are very different from the "Drude" profile. Fortunately, the 
physics of carbon grain formation and evolution apparently 
does not allow the evolution towards such extreme shapes. 
Besides, the grain surface is certainly very irregular and may 
not be subjected to the same theoretical treatment as per- 
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Figure 16. Effect on the UV bump of hydrogenation (H/C) of the 
faces of graphene stacks, red: clean faces; green, blue and black: 
H/C(per face)=50, 75 and 100%, respectively. The corresponding 
widths are 1, 1 and 0.65 iim~^] and the peak wavenumbers: 4.63, 
4.61 and 4.84 /im-^. 



fectly smooth ellipsoids. Finally, even though the term "sur- 
face mode" is usually applied to the Frohlich resonance, it is 
clear that the whole grain volume is involved, so the surface 
shape may not be so important, after all. In order to sub- 
stantiate these conjectures, we also computed Q/a for values 
of L symmetrically bracketing 1/3, from 0.25 to 0.4, corre- 
sponding to ratios of principal axes ranging between 1 and 
about 2. It was found that, although the peak wavenumber 
of the feature increases from 4.4 to 4.82 /im""*^, the average 
feature is still narrow (1 /im~^) and does not shift much 
away from 4.68 /im~^. This is due to the symmetry of the 
variations, a consequence of the implicit assumption that 
the fluctuations of L around 1/3 are limited and symmetric. 



7.5 Turbostratic disorder 

Inter- atomic forces naturally tend to maintain adjacent 
graphene layers in the preferred stacking schemes, ABA 
or ABC. Recently, however, slightly disoriented graphene 
multi-layers have attracted attention in the laboratory (see 
Kim et al. 2012; Trambly de Laissardiere et al. 2010). We 
also considered the impacts of this type of disorder upon 
the dielectric function. For this purpose, we performed pre- 
liminary calculations, using a simplified tight-binding the- 
ory (assuming constant electric dipole moment). The high 
energy wing of the tt resonance was found to be quite in- 
sensitive to disorientation, and so is the UV bump. But, 
interestingly, the resonance peak is split into 2 peaks, whose 
separation increases with the angular deviation (^ 1 eV for 
10 deg between adjacent layers). This effect could contribute 
to the flattening of the resonance peak of terrestrial and lab- 
oratory graphite (see Taft and Philipp 1965). 



7.6 Hydrogenation 



Duley and Seahra ()1998p considered in detail the effects of 
edge hydrogenation on various regular aromatic structures. 
They noted some positive correlations between 7 and hydro- 
genation (their Fig. 11). For partial hydrogenation, there are 
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also some weak positive correlations between uq and hydro- 
genation. 

In our treatment, we use rectangular graphene sheet 
with armchair and zigzag edges. The size of a sample with 
4096 X 4096 carbon atoms is about 8700 x 5000 ^ 0.44/im^ 
This is large enough that its electronic properties are largely 
insensitive to the (open or periodic) boundary conditions 
(see Yuan et al. 2010). We assume, here, that this is also 
true for the flakes we consider in the dust model of Fig. 11 
(^ 100 X 100 A^ in size). It must be kept in mind, of course, 
that edge hydrogenation of large PAHs was shown to affect 
notably tt — tt* transitions (Duley and Seahra 1998; Malloci 
et al. 2008). However, the maximum size of PAHs considered 
in these studies was smaller than 7x7 rings, which is way 
smaller than our assumed model dust. 

Internal random hydrogenation considered in the 
present work is totally different from edge hydrogenation in 
that it introduces sp^ hybridization, but it was also found 
to increase both u (slightly) and 7 (see Fig. [M]) . 

Another possible "defect" is the hydrogenation of the 
outer faces of the edge layers of stacks of graphenes. The ef- 
fect on the bump is illustrated in Fig. [16] for 3 degrees of hy- 
drogenation: H/C=50, 75 and 100 %. It appears that partial 
hydrogenation only broadens slightly the bump. Complete 
hydrogenation, however, seems incompatible with observa- 
tions. This may mean that hydrogenation in space is only 
partial. But this may also be due to the fact that our tight- 
binding model for hydrogenated graphene only includes the 
TT band, and is therefore not accurate for large concentration 
of hydrogenation, in which there is a strong overlap between 
the a and tt bands. A tight-binding model including both a 
and TV bands for hydrogenated graphene is currently under 
development. 

By contrast, the absence of the a bands in our present 
treatment is not expected to affect the behavior of graphene 
upon the introduction of the other two types of disorder 
considered in this work (random on-site potential, random 
hopping constant); for the latter do not add any extra terms 
to the Hamiltonian. 



7.7 Ionization 

Ionization of the grains is likely to occur in the ISM and 
should therefore be considered here. We found that an elec- 
tronic charge in the order of 0.1 electron per C atom will 
only shift the Fermi level by 0.2 eV, leaving the spectrum 
at higher energies unchanged. In fact, taking into account 
that the tt resonance is associated to inter-band transitions 
from the Van Hove singularity of the valence band to the 
Van Hove singularity of the conduction band, both of them 
being separated by an energy window larger than 4eV, weak 
amounts of doping will not modify the optical spectrum 
enough to affect the main features of the resonant peak of 
the TT mode, (see Yuan et al. 2011a). 

On the other hand, Cecchi-Pestelhni et al. ()2008p 
showed that, for several PAH's, the spectra of neutrals, an- 
ions, cations and dications were different between 6 and 10 
eV. This effect may be due to interaction between the tv 
and a resonances. Addressing this issue will become possi- 
ble when both resonances are included in the tight-binding 
model. 



7.8 Work in progress 

Clearly, a valuable addition to this work would be the in- 
clusion of the a resonance in the calculation of the DOS 
and the optical conductivity. Work is in progress to fill this 
gap. To include both a and tv bands one needs to consider 
four orbitals (s, px, py and pz) on each carbon atom and 
one orbital (s) on each hydrogen atom. The dimension of 
the hopping matrix between two nearest carbon atoms be- 
comes 8x8, which makes the computation more expensive 
as compared to the single tt band model which only includes 
a 2 X 2 hopping matrix. Each 8x8 hopping matrix between 
two neighbours has to be recalculated according to the co- 
ordinates of each carbon atom. To the authors' knowledge 
the full tight-binding model including both a and tv bands 
for graphene is not available in the literature yet. 

At the same time, we are investigating whether the ac- 
curacy of the present tight-binding calculation can be im- 
proved, for instance, by using larger samples to have more 
random complex coefficients in the initial wave function, or 
using more time-steps to increase the energy resolution in 
the spectrum. The fundamental operation in our numerical 
methods is the time-evolution of the wave function by using 
the Chebyshev polynomial method, which has already the 
same accuracy as the precision of the simulation machine. 



8 CONCLUSION 

The present model carrier consists of nano-sized "bricks", 
i.e. stacks of a few finite sheets of graphene. The physics 
of carbon is now known in such detail that it has become 
possible to compute the dielectric functions of such particles 
from first principles, with no ad hoc assumptions, nor tailor- 
ing of any sort. In its clean and ordered form, this structure 
carries a UV bump at 4.65 /im"\ with FWHM 0.85 /im"\ 
very near the lower limit set by observations. This width 
is essentially determined by the high-energy slope of 62 as 
a function of frequency, which, in turn, depends heavily on 
the variation of the electric dipole moment with frequency. 
Unlike theories which assume for simplicity that the dipole 
moment is constant, the present treatment explicitely in- 
cludes its variations. Because the dipole moment decreases 
notably as frequency increases, the slope of £2 becomes no- 
tably steeper, and the bump width distinctly narrower. 

Most natural or artificial carbonaceous materials are 
neither pure nor perfectly regular in structure. Impure or 
disordered structures can also be simulated theoretically by 
allowing the natural energetic parameters of pure graphene 
to vary according to a statistical distribution, thus simu- 
lating different types of structural, chemical and orienta- 
tion disorder. We have found that three types of struc- 
tural/chemical disorder affect notably the high-energy slope 
and, hence, broaden the UV bump. Moderate degrees of dis- 
order are enough to cover the observed range of widths while 
the peak position remains constant to ^ 1 %. 

The peak position is more sensitive to the type of dis- 
order than to its degree. For any of these types of disorder, 
the general trend is for the width to increase with the degree 
of disorder. Orientation disorder of bricks randomly glued 
together in one dust particle, brings about more broaden- 
ing and peak shift. This dust model, too, is sensitive to the 
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structural and chemical disorders considered in this paper, 
with similar consequences on position and breadth of the 
peak. When all cases are lumped together in the same dia- 
gram, no obvious correlation emerges between bump width 
and position. 

In essence, our theoretical procedure boils down to re- 
placing sp^ by sp^ or sp^ bonds. No wonder, then, that this 
leads, in particular, to a depression of the tt resonance, which 
is observed in various amorphous carbons and hydrocarbons 
or coals of various ages. Together with the general trend of 
increasing width and decreasing strength of the resonance, 
the results that we present here can help to draw parallels, 
and identify common features, between the different pro- 
posed bump-carrier models. 
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